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[^^ ; ABSTRACT 

Aims. Observational data on companion statistics around young stellar systems is needed to flesh out the formation 
pathways for extrasolar planets and brown dwarfs. Aperture masking is a new technique that is able to address an 
important part of this discovery space. 

Methods. We observed the two debris disk systems HD 92945 and HD 141569 with sparse aperture masking (SAM), a 
rG • new mode offered on the NaCo instrument at the VLT. A search for faint companions was performed using a detection 

Mh' strategy based on the analysis of closure phases recovered from interferograms recorded on the Conica camera. 

Results. Our results demonstrate that SAM is a very competitive mode in the field of companion detection. We obtained 

5ct high-contrast detection limits at X/D of 2.5 x 10"^ (AL' = 6.5) for HD 92945 and 4.6 x 10"^ (AL' = 5.8) for 

HD 141569. According to brown dwarf evolutionary models, our data impose an upper mass boundary for any companion 

for the two stars to, respectively, 18 and 22 Jupiter masses at minimum separations of 1.5 and 7 AU. The detection 

limits is mostly independent of angular separation, until reaching the diffraction limit of the telescope. 

Conclusions. We have placed upper limits on the existence of companions to our target systems that fall close to the 

planetary mass regime. This demonstrates the potential for SAM mode to contribute to studies of faint companions. 

^ ' We furthermore show that the final dynamic range obtained is directly proportional to the error on the closure phase 

\^ ' measurement. At the present performance levels of 0.28 degree closure phase error, SAM is among the most competitive 

fSJ , techniques for recovering companions at scales of one to several times the diffraction limit of the telescope. Further 

^+ . improvements to the detection threshold can be expected with more accurate phase calibration. 

Key words. Instrumentation: high angular resolution - Stars: planetary systems - Stars: individual: HD 92945 - Stars: 
t^ ■ individual: HD 141569 
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'^ 1. Introduction 8 and 15 Au) using direct deep imaging validates the link 

J> , bet ween structures in debris dis ks and the presence of plan- 

V~j ■ Since the first iniaging of a debris disk around /3 Pictoris ets (jLagranee et al.l[2009ll2010t) . Besides this planet, only a 

rS I (jSmith &: Terrild ll984D . and the remarkable successes of couple of planetary- mass objects have been detected around 

H the radial velocity technique for exoplanet discovery her- stars with a debris disk. In fact, most debris disks have been 

■ - - ■ aided nearl y a decade later by a companion to solar-like imaged around young bright A- or F-type stars, which are 

51 Pegasi (jMavor fc Quelod [l995D . the understanding of either rapid rotators and/or active stars. The search for 

planetary systems formation and evolution has become one companions with radial velocity techniques around these 

of astronomy's greatest challenges. Today, at least 22 de- stars is therefore extremely difficult. Key breakthroughs 

bris disks have been imaged at optical, infrared, or sub- have recently been made with high-contrast imaging: a < 3- 

milhmeter wavelengths. Such disks are thought to be the Mjup planetary companion was detected in the outskirts of 

cradle where planetary systems have recently formed. The Foma lhaut's debris disk (119 AU from the star: lKalas et all 

observed dust component is not primordial but instead pro- 120080 . while four planetary companions of 7, 10, 10, and 

duced by collisions among larger rocky bodies. This makes 10 Mjup were im aged at 68, 38, 24, and 14 AU, respectively, 

debris disks obvious places to search for planets, and many from HR 8799 (jMarois et all I2008L |2010[ ). It is not clear 

show possible indirect signs of their presence, such as dust whether these objects could have formed through core ac- 

structures or disk asymmetries. cretion - and not from gravitational instabilities - the fore- 

The recent detection and confirmation of a - 9 Mjup most theory of planet formation, although /3 Pic b is a se- 

planet within the disk of /3 Pic (semi-major axis between 
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Fig. 1. Left panel: Diffraction pattern observed on star HD 135549. This image is a single frame of 0.4 ms integration 
time, using the 7- hole pupil mask. Background noise is ~12 ADU, and pic brightness 800 ADU. Right panel: Fourier 
transform of the diffraction pattern. Each peak corresponds to a baseline pair. 



rious candidate for this formation scenario ([Lagrange et al.l 
[2OIOI) . 

Most companions detect ed with h igh-contrast iniag- 
ing, such as AB Pic b (Chauvin et al.l i2005') or 2M1207 b 
(|Chauvin et al.ll20M ). are believed to be either too far or 
too massive to have formed through core accretion. An 
alternate pathway gaining widespread interest is a stellar 
formation mechanism, i.e., formation during the fragmen- 
tation and collapse of a molecular cloud. There are many 
unanswered questions raised by these models. Do these two 
processes for forming planets and brown dwarfs/low-mass 
stars operate exclusively (does one inhibit the other)? Can 
core accretion from a disk operate when a brown-dwarf 
companion exists at large separation? In contrast, can disks 
and planets exist around tight binary systems, composed of 
the primary star and a very close low-mass star or brown 
dwarf? How do massive companions impact the dynamic 
evolution of planets and disks? 

To answer these questions, it is necessary to probe the 
innermost regions of debris disks when searching for sub- 
stellar companions. For the closest debris disk, direct imag- 
ing proved to be a successful technique, as shown for /3 Pic 
or Fomalhaut. But most of the young associations (as well 
as stellar forming regions) lie at more than a hundred par- 
sees. At such a distance, 10 AU is equivalent to 100 mas, less 
than twice the diffraction limit of an 8-meter telescope in 
the near-infrared. This observational domain is not filled by 
differential imaging techniques, whose inner working angle 
is typically of few resolution elements. It is neither filled 
by long baseline interferometry: the field of view of a single 
mode interferometer is too small, equal to the resolution 
element of one of its telescopes. The effectiveness of b oth 
techniques is reviewed in detail in lAbsil fc Mawetl(|2010t ). In 
terms of angular separation, aperture masking is the miss- 
ing link in the field of high dynamic observations. It gives 
a relatively high dynamic range (5 to 7 mag) at a specific 
angular separation (sa 100 mas) which is of greatest interest 
for planetary formation. 

In this paper, we present three stars observed with an 
aperture mask at the VLT. In the next section, the instru- 
ment is described, along with the observational procedure 
and data reduction methods. In section [31 we present a 



clearly resolved detection of the binary star HD 135549 for 
illustrative and pedagogical purposes. Section [3| reports on 
the non detection contrast limits obtained with SAM for 
the two debris disks HD 92945 and HD 141569. We present 
these limits as a function of flux ratio and infer the mass 
limit derived from current models. 



2. Observation &. data reduction 

2.1. SAM at NaCo 

An aperture mask works by transforming the telescope into 
a Fizeau interferometer. The point spread function is a 
complex superposition of fringes at given frequencies. In 
specific cases, pupil masking can outperform more tradi- 
tional differential imaging for a number of reasons. First, 
the masks are designed to have nonredundant array config- 
urations that permit phase deconvolution; slowly moving 
optical aberrations not corrected by the AO can be accu- 
ratly calibrated. Second, the mask primarily rejects base- 
lines with low spatial frequency and passes proportionately 
far more baselines with higher X/B resolution than does 
an orthodox fully filled pupil. Third, high-fidelity recov- 
ery of phase information allows "super resolution" , with a 
marginal loss of dynamic range up to A/2D. The principal 
drawback is a loss in throughput so that photon and de- 
tector noise can affect the signal-to-noise ratio even where 
targets are reasonably bright for the AO system to close 
the loop. In common with all synthesis imaging by inter- 
ferometry, the effective field-of-view of SAM is determined 
by the shortest baseline so that the technique is not com- 
petitive at separations larger that are greater than several 
times the formal diffraction limit. 

VLT's SAM mode has grown from earlier seeing- limited 
experiments on the Keck telescope (Tuthill et al. 2000), and 
more recently has ga ined widespread use on adaptive optic s 
equipped telescopes ([Tuthill et al.ll2006t[Llovd et al.ll2006D . 
Commissioning of SAM at NaCo has required important 
software development from ESO at Paranal. It prompted 
the development of datacube storage (offered now regularly 
on NaCo as "cube mode") and of pupil-tracking observa- 
tion (the Nasmyth rotator tracks the telescope elevation). 
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SAM is currently under investigation as an additional mode 
for the VISIR instrument, operating in the mid- infrared. 
It is under inform al study for the SPHERE instrument 
(jBeuzit et al.ll2008D . 



2.2. Observing procedure 

Observations were taken in the L prime band (A = 3.80 ± 
0.31//m) using the "7 holes" mask and the IR wavefront sen- 
sor (WFS). The diffraction pattern is plotted in Fig.[T] The 
physical properties of the mask ~ as well as the three others 
- are described in the NaCo manuaQ- The complete obser- 
vational sequence for an object typically took two hours. 
Table [T] contains observing details such as integration times, 
datacube size, and number of acquisitions. 

Each science target is associated to a PSF reference star 
situated within one degree. A single observing block (OB) is 
used that encodes replicated observations for both calibra- 
tor and target. Rapid star acquisition - which is important 
for ensuring good calibration - can be augmented by using 
a simple offset of the telescope without having to preset the 
telescope and perform a full (instrumental) re-acquisition. 
In this mode, which has been christened "star hopping" , the 
adaptive optics loop is opened, while the template orders 
a dither to bring a different star into the AO field selector. 
The AO is then closed manually by the operator without in- 
curring the time penalty for re-optimization. Star hopping 
therefore only works on objects of comparable brightness 
in the wavefront sensor. 

As the observation progresses, the repetition of the tem- 
plate collects eight datacubes of multiple frames (typically 
a hundred), each at a different dither position on the de- 
tector. For these observations, the detector was windowed 
to 512 X 514 pixels. 

A single snapshot observation yields Fourier coverage 
of 21 spatial frequencies. However, the pupil-tracking mode 
results in sky rotation on the detector as the parallactic 
angle changes. This variation with time results sweeps the 
baselines into circular Fourier tracks and permits rotational 
aperture synthesis techniques to assist with the filling of the 
spatial frequency UV-plane plane, as illustrated in Fig. [21 

2.3. Data reduction 

To the author's knowledge, two data reduction software li- 
braries are presently availablqj to reduce SAM data. One 
pipeline has been developed by Sydney University, Cornell 
University, and Caltech from 2004 onwards, based on an 
earlier pipeline from Berkeley. It has already been used 
for several papers arising predominan tly from Keck aper- 
ture masking data (JTuthill et al.l 120001 ). This reduction al- 
gorithm is based on fast Fourier transform (FFT). The 
data presented in this paper have been reduced by the 
sparse aperture mode pipehne (SAMP) developed at the 
Observatoire de Paris. This software is similar to that from 
Sydney university, and numerous tests and cross-checks 
have produced similar results. 

For both pipelines, data reduction follows a similar 
path: 

^ http://www.eso.org/sci/facilities/paranal/instruments/naco 
^ upon request to developers 
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Fig. 2. Frequency coverage 
HD 92945, and HD 141569. 
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Fig. 3. Phase as a function of the spatial frequency. The 
different plots are models of a binary system with differ- 
ent flux ratio and angular separation. No noise has been 
added, so the curves represent the perfect case. From top 
to bottom, the contrast of the primary to secondary de- 
crease. Panels (a-b-c): flux ratio of 1%. Panels (d-e-f): flux 
ratio of 0.1%. Panels (g-h-i): flux ratio of 0.01%. From left 
to right, the angular separation gets smaller. Panels (a-d- 
g): 200 mas. Panels (b-e-h): 100 mas. Panels (c-f-i): 50 mas. 
In the L' band, the cut-off frequency of an eight meter tele- 
scope is 0.01 mas^^. 
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Table 1. Observation log 



Target 


Date 


Start time (UT) 


End time (UT) 


DIT (ms) 


NDIT 


Dithers 


Acquisitions 


HD 135549 
HD 135344b 


2010-03-18 
2010-03-18 


05:43:27 
05:32:31 


10:13:37 
10:04:16 


0.2 
0.2 


200 
200 


8 
8 


15 
15 


HD 92945 
HD 92933 


2010-04-07 
2010-04-07 


02:06:41 
02:14:49 


03:01:27 
02:52:55 


0.4 
0.4 


100 
100 


8 
8 


4 
3 


HD 141569 
BD-03 3826 


2010-04-07 
2010-04-07 


07:16:35 
07:34:55 


08:33:03 
08:24:37 


1.2 
1.2 


35 
35 


8 
8 


4 
3 
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Fig. 4. Theoretical (Monte-Carlo) 1 a detection limits of 
the 7-hole mask as a function of angular separation and 
error on closure phases. The ripples on the detection limit 
curves are not due to uncertainty in the Monte-Carlo sim- 
ulation, but to the nonhomogeneous coverage of the UV 
frequency plane. 



— Pre-reduction: the data are (i) fiatfielded, (ii) bad- 
pixel corrected, and (iii) background-subtracted. 

— Complex visibility: Complex visibilities are obtained 
from FFT of the image (Sydney pipeline), or by direct 
fringe fitting (SAMPfl. 

— Bispectrum: the bispectrum is obtained using multi- 
plication of the complex visibilities corresponding to 
a triangle of holes in the mask. Bispectrum is aver- 
aged over all the frames of each set of eight datacubes. 
Closure phases are obtained using the argument of the 
bispectrum. 

— Calibration: Closure phases obtained on the calibrator 
are fitted with a polynomial function, and subtracted 
from the closure phases on the science object. 

— Detection: The closure phases are fitted with a model 
of binary star (see section [3]). A x^ map is computed, 
giving the likelihood of the presence of the companion 
for any given separation and position angle. 

Atmospheric dispersion is not accounted for in our 
pipeline, but we are aware that in certain conditions (pres- 
ence of emission lines in the spectrum) , it can mimic asym- 
metry along the elevation axis. Such false detection can be 
avoided by observation at different parallactic angles. 



The effect of fringe blurring due to the large bandwith 
of L' filter can technically be accounted for in the fitting 
process. It would require including a loss of fringe coher- 
ence on the edge of the diffraction pattern. However, at 
this stage, the SAMP pipeline does not include this correc- 
tion. It is, however, not seen as critical with the seven-hole 
mask, since the diffraction pattern of the 1.2 meter holes are 
significantly smaller than the coherence length (see Fig.[T]). 
It becomes an acute problem with the small holes of the 
18-hole mask, which then justifies using a narrow band fil- 
ter. 

The complex visibilities consist of amplitudes and 
phases, but for the purposes of binary detection, we typi- 
cally rely on the phases only (of the bispectrum) . This is be- 
cause of the dependence of absolute visibilities on the adap- 
tive optics correction. They are thus harder to calibrate, 
and are subject to different biases between the science and 
the calibrator object. However, by default, the pipelines 
also process the power spectrum and estimate fringe visi- 
bilities. Performance of SAM's visibility estimation will be 
the subject of a forthcoming paper. 



3. Test case: the HD 135549 binary system 

3.1. Principles of faint companion detection witli SAM 

Aperture masking data yields Fourier phases and ampli- 
tudes. They are functionally equivalent to those recov- 
ered in long-baseline interferometry, and therefore the same 
techniques and software can be applied. Image reconstruc- 
tion is a possibility, but to fulfill the difficult requirements 
of faint companion detection, direct fitting in the Fourier 
domain is essential. 

Modeling the presence of a two point-like sources in 
the Fourier domain is trivial. The complex visibility can 
be written as: 



V{u, v) — 1 ~ r {1 — exp [i2Tr{—au + Sv)/X]} 



(1) 



^ Fringe fitting - in a least squares sense - is performed 
through matrix inversion thanks to singular value decomposi- 
tion 



where u and v are orthogonal spatial frequency vectors (re- 
spectively to the east and to the north), a and S are the 
right ascension and declination in radians, and r is the fiux 
ratio between the primary and secondary. Figure [3]presents 
the interferometer fringe phase for a grid of nine models 
with three different angular separations and three different 
flux ratios. A direct relationship between the phases and 
the basic properties of the binary system can be noted. The 
phases exhibit a sinusoidal curve across the Fourier plane, 
with the amplitude proportional to the contrast ratio while 
the period is proportional to the separation. Closure phases 
are derived from the baseline phases as a simple summa- 
tion of the three phases corresponding to a closed triangle 
of holes. 
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RA [mas] 

Fig. 5. x^ surface as a function of a and 5 obtained by 
fixing r to the best-fit value in the model of a binary stated 
in Equation [TJ A clear minimum indicates the position of 
the stellar companion. 



Companion searches are performed by fitting the vari- 
ous binary parameters (r, a, and (5) to the measured clo- 
sure phases, so the detection limit obtained from any given 
dataset is directly related to the signal-to-noise ratio of the 
closure phase. A Monte-Carlo simulation was performed to 
test the theoretical detection limit, as a function of the 
angular separation and closure phase measurement error 
bars. The result is presented in Fig. 21 The detection limit 
is almost constant for any separation larger than A/D. The 
Monte-Carlo simulation shows that the dynamic range can 
be approximated as (as a function of a{CP) in degrees) 



detection (la) = 2.5 x 10"^ x a{CP) . 



(2) 



Detections below A/D can still be significant, but the con- 
trast achieved decreases rapidly with the separation. At half 
A/D, limits are worse by a factor of two, while at 0.25 A/D 
the degradation in performance is a factor 15. For separa- 
tions beyond the formal diffraction limit, any fitting process 
becomes highly model dependent and often strong degen- 
eracies between fitted parameters can appear. A typical ex- 
ample at small separations is the covariance between the 
flux ratio and separation: it becomes difficult to measure 
these two parameters independently. 

3.2. Detection of the stellar companion to HD 135549 

HD 135549 was initially observed for the purpose of a PSF 
reference star, but it turned out to be a binary system with 
a contrast ratio around 5% and a separation close to 180 
mas. We present the data in this paper for the sake of in- 
troducing the technique of detecting faint companions with 
SAM illustrated by strong, clear systematic signals. 

The free parameters of Eq. ([Ij are the binary coordi- 
nates a and S, along with the flux ratio r. A map of the 
reduced x^ is derived as a function of a and S (Fig. [5]). It 
indicates the presence of a companion with a strong likeli- 
hood at a separation of « 180 mas to the southeast of the 
star. The best-fit model is compared to the observed data 
in the two panels of Fig. [Sj The lefthand panel shows the 
observed versus model-predicted closure phase data. The 
final reduced x^ is equal to 1.85 while the average residual 
is 0.28 degree over a total of 35 x 15 = 525 closure phase 
measurements. According to Eq. ([2]), this dataset therefore 



has the capability of a detection limit of 7.5 x 10^^ (1 a de- 
tection) . The righthand panel shows phases recovered from 
the closure phases. This is done by inverting the phase to 
closure phase matrix. The matrix is not fully invertible; for 
example, tip-tilt is not constrained by the closure phase. 
Therefore, unknown parameters of the inversion were taken 
to match the model. The best fit model is over-plotted as 
a solid line. 

The first 15 lines of Table [5] results from the processing 
of each SCI-CAL acquisition pair separately. An indepen- 
dent model is fitted at each epoch, and the detection can be 
seen to be rotating in the detector frame (but fixed on the 
sky). A strong indication that our detection is indeed real 
lies in the rotation in the detector frame following the par- 
allactic angle evolution. Because the telescope is azimuthal, 
it is unlikely that a calibration error (including optical aber- 
rations) rotate in such a motion. This validation could not 
have been made with an equatorial telescope. 

The last line of Table [5] gives the best- fitting binary pa- 
rameters. Separation and dynamic ranges are obtained with 
high statistical precision. The precision of the closure phase 
measurement over time is highlighted in Fig. [71 in which 
each panel depicts a different closure phase. The solid line 
corresponding to the model can be seen to almost perfectly 
match the evolution of the closure phase as a function of 
time. 

3.3. On the limitation of the dynamic range 

As stated by Eq. [21 the dynamic range is a pure function 
of the error on the phase. The binary system was observed 
during good atmospheric conditions (0.7" seeing and a co- 
herent time of 8 ms) with a relatively fast frame time of 
200 ms. Phases recovered from the interferograms show a 
median frame-to-frame error of 17 degrees, equivalent to a 
180-nm piston jitter at our observing wavelength. On the 
other hand, the frame-to-frame closure phase error is only 
2.8 degrees. This is a clear indication that most of the 180- 
nm rms error is due to the residual of atmospheric pertur- 
bation not corrected by the AO system. 

Assuming that these residuals were caused by photon 
noise or readout noise, the remaining 2.8 degrees frame- 
to-frame error should be uncorrelated. In this case, the fi- 
nal error on a whole acquisition sequence (8 dither of 200 
frames) should improve as the square root of the number of 
frames giving an expected error of 0.07 degrees. For such a 
scenario, the dynamic range obtained in the final detection 
should be 1.7 x 10~* (according to Eq. ^). Unfortunately, 
the observed standard deviation for both calibrator and the 
science target closure phases are significantly larger (0.25 
and 0.33 degrees respectively; obtained by comparing dif- 
ferent telescope dithers). These larger errors point to an 
additional systematic term generating biases in the phase 
closure which do not average out over long data runs. In 
order to extract the very best contrast ratio in our faint 
companion detection strategy, these closure phase biases 
need to be understood and calibrated. 

Within our present software pipeline, closure phase cal- 
ibration is performed by subtraction of the PSF reference 
star closure phases which resulted in an rms closure phase 
error observed of 0.24 degrees (when comparing different 
dither positions). Compared to the theoretical limit of the 
high-frequency noise (v^ x 0.07 w 0.1 degrees), it is appar- 
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Fig. 6. Lefthand panel: Measured closure phase as a function of model closure phase assuming the best-fit binary 
parameters. The rms error between the model and the measurement is 0.28 degrees. Righthand panel: phase obtained by 
least square inversion of the closure phase to phase matrix, plotted as a function of the projected baseline (perpendicular 
to the binary orientation; hence the sinusoidal curve). 




Hour (UT) 

Fig. 7. Closure phase as a function of time for the 10 largest triangles. The triangles are made of the holes (a) 1-5-7, 
(b) 1-2-6, (c) 1-6-7, (d) 4-5-6, (e) 1-2-7, (f) 1-5-6, (g) 4-5-7, (h) 1-4-7, (i) 2-5-7 and (j) 2-5-6; numbered according to the 
rows of Table 4-19 in the NaCo manual (available online at eso.org). Calibrator data are represented as circles, while 
binary star data are squares. The dashed line is the polynomial fit of the calibrator, and the solid line is the closure phase 
predicted by the best fitting binary system model. The same fit quality is obtained on all 35 closure phases. Overall, the 
average rms error of the fit is 0.28 degrees. 
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Table 2. Detection as a function of time using individual pairs of source-calibrator files. 



UT 


Position Angle 


Separation 


Position Angle 


Flux 




on detector (deg) 


(mas) 


on sky (deg) 


(% of central object) 


05:43:27 


-153.38 


179.32 ± 0.38 


148.24 ± 0.13 


5.62 ±0.07 


06:03:35 


-148.63 


179.33 ± 0.37 


147.96 ±0.11 


5.51 ±0.07 


06:21:57 


-143.54 


181.52 ±0.26 


148.25 ± 0.09 


5.69 ± 0.04 


06:40:14 


-137.46 


179.55 ± 0.58 


148.64 ±0.21 


5.60 ± 0.07 


06:58:02 


-130.17 


179.46 ± 0.37 


148.28 ±0.10 


5.62 ±0.06 


07:16:34 


-120.66 


179.24 ± 0.30 


149.15 ±0.12 


5.60 ± 0.05 


07:35:08 


-108.48 


179.42 ± 0.29 


148.41 ±0.10 


5.57 ±0.06 


07:53:06 


-93.93 


179.34 ± 0.46 


148.59 ±0.18 


5.45 ± 0.08 


08:11:23 


-77.02 


180.50 ±0.91 


148.10 ±0.21 


5.43 ± 0.08 


08:29:18 


-60.32 


180.80 ± 0.50 


148.41 ±0.19 


5.47 ±0.07 


08:55:48 


-39.29 


181.44 ±0.28 


148.72 ±0.10 


5.58 ±0.04 


09:13:54 


-28.29 


180.80 ± 0.30 


148.50 ±0.12 


5.57 ±0.05 


09:37:18 


-17.41 


181.60 ±0.39 


148.56 ±0.11 


5.27 ±0.07 


09:55:46 


-10.75 


180.40 ± 0.54 


148.07 ±0.16 


5.33 ±0.07 


10:13:37 


-5.44 


179.73 ± 0.34 


147.61 ±0.13 


5.49 ± 0.04 


Full Fit 


all 


180.02 ±0.11 


148.36 ± 0.04 


5.53 ± 0.02 



ent that systematic biases in the closure phase remain the 
dominant limit to the dynamic range. When fitting the fi- 
nal data with a model of a binary system, the residual error 
on the closure phases is 0.28 degrees. 

This statistical analysis shows that the high-frequency 
noise - assumed to be photon and readout noise - is not 
limiting the dynamic range. On this specific dataset, one 
magnitude in contrast ratio could be gained if only it was 
possible to have a perfect calibration of the closure phase 
systematics. A higher dynamic range could then be achieve 
by repetition of the same observation. The reason for these 
systematics errors is still under debate. Between several 
possibilities, we can mention i) optical aberrations before 
the pupil mask, ii) correlated detector read out noise, iii) 
spatially correlated phase error due to the modal correc- 
tion of the adaptive optic system. The residuals could be 
diminished partially by observing several calibrators. But 
not to a full extent: the noise would average only on the 
calibrator dataset. 



4. The two debris disks HD 92945 and HD 141569 

4.1. The targets 

4.1.1. HD92945 

HP 92945 (V 419Hya) is a K1.5V star located at 22 pc. 
(jSong et al | [2004) estimated an age of 100 Myr using Can 
lines. I f we suppose HD 92945 to b e part of the AB Dor 
group (jLopez-Santiago et al.l 120061 ) . this age agrees with 
the estimation of ~ 50-120 Myr for the whole group 
(jZuckerman et all2004l : iLuhman et al.ll2005[ ). Observations 
with the Advanced Camera for Surveys (ACS) on the 
Hubble Space Telescope have revealed it is surrounded by 
a disk ^ 30° from edge-on with a O'.'7-thick bright ring 
at 57 AU from the star. The disk has a diffuse componen t 
detected between 55 and 170 AU (jGolimowski et al.ll2007t) . 
The inner gap may well be the res ult of resonances with 
a ste llar or planetary companion (jArtvmowicz fc Lubowl 
119941) . 



Table 3. Best-fit detection 



Target 


Degrees of 
freedom 


Single star 


Min x^ 
Binary model 


HD 135549 


522 


248318 


908 


HD 92945 


137 


282 


234 


HD 141569 


137 


265 


223 



4.1.2. HD 141569 



HD 141569 is a 5 ± 3-Myr-old (jMerin et al.ll2004D pre-main 
sequence B9.5V star in a triple system. It is located 99-pc 
away and the central star features an ele gant debris disk 
that was first imaged in the near-infrared (Augcreau et ajj 
1999) and in visible scattered light dClampin et al. 2003J ). 
The two other components of the syste m form an M2-M4 bi- 
nary, ~ 7'.'5 away from the central star. iReche et al.l (|2009tl 
give an overview of the rich disk structure. It is composed 
of two annuli with peak luminosity in scattered light at 
~ 200 and ^ 325 AU from the star. These two rings have az- 
imuthal brightness asymmetries (factors 2 to 3). The outer 
ring shows a tightly wound spiral structure. The rings are 
separated by an apparent gap wider than the rings. An ex- 
tended diffuse emission associated with a faint spiral arm is 
present in the northeast side of the disk, and is detected to 
more than 600 AU. Another spiral arm is possibly observed, 
pointing toward the binary companions. Finally, the disk 
brightness sharply decreases within 200 AU, suggesting a 
strong dust depletion in the innermost regions of the disk. 

4.2. 5a detection limits 

Establishing a high-contrast detection limit follows a differ- 
ent approach than establishing a detection. The most im- 
portant step is indeed to show that there is no significant 
detection in our data. A practical approach is to analyze 
the x^ maps (Fig. [5]). If the minimum value of the x^ does 
not significantly differ from the maximum x^ i it means that 
the binary model does not significantly improve the fit to 
the data. The x^ values at best fit are given in Table [31 
The detection for the binary system HD 135549 discussed 
in the previous section is extremely clear with a normalized 
X^ of 1.85. This dramatically improves the fit compared to 
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Fig. 8. x^ as a function of a and S obtained from the best fit of r in the model of a binary stated in Equation [TJ Left 
panel: HD 92945. Right panel: HD 141569. The lower limit of the color scale (139) is the number of degree of freedom of 
the fit. 
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Fig. 9. x^ maps for target HD 92945 (left panel) and HD 141569 (right panel). The solid curves correspond to the 
isocontours at min(x^)+l; +4; and +9 (eg. one, two, and three sigmas contours). Multiple minimums are present (false 
detections). The 5a detection limits are, averaged over the separation, 2.5 x 10~^ and 4.6 x 10"'^ (dashed curve). The 
right vertical axis was obtained assuming a linear fit of DUSTY evolutionarily models. 



the null hypothesis (no binary x^ /"-freedom = 423). This is 
not the case for the two debris disk systems. At best fit, 
the binary model improves the x^ by a factor of 1.24 for 
HD 92945 and 1.18 for HD 141569. 

Although the x^ maps indicate a nondetection, robust 
statistical limits based on these data remain to be derived. 
To accomplish this, we plot the x^ in Fig. [HI minimized 
over position angle, as a function of separation (y/a^ + S^) 
and contrast (r). The red solid contours correspond to 1, 
2, and 3 a detection limits. Several minima exist in the 
maps including two 3 a false detection: at 440 mas (contrast 
0.16%) and at 160 mas (contrast 0.2%). It is noteworthy to 
point out that the absence of a companion (equivalent to a 
binary system with a separation of mas) is within 5 cr of 
the false detections. This 5 a limit corresponds to the red 
dashed curves. We define this curve as our 5 a detection 
limits, hence excluding any detection on the two dataset. 
This approach is coherent with the generally recognized 5a 
criterion to validate a detection. 

According to Fig. ^ the average 5 a detection limits of 
HD 92945 and HD 141569 over the separation range 100- 
500 mas are respectively 2.5 x 10~^ and 4.6 x 10~^, that 
is, AL' of 6.5 and 5.8. However, the nondetection limit is 
not completely uniform over the field of view and can have 



values that are nearly a factor of two worse. We can ob- 
serve the detection limit worsening by a factor 2 between 
X/D and 0.5 X/D, with a precipitous drop for still lower 
spatial separations. These results also confirm the Monte- 
Carlo simulations (presented in Fig.U]), assuming error bars 
on the closure phase of 0.2 and 0.4 degrees. 



4.3. Jupiter-mass limits 

To convert contrast limits into upper bounds for the mass of 
any possible companion, we must know the absolute mag- 
nitude corresponding to the detection limit in each system. 
We derived the L' band magnitudes from the K band mag- 
nitudes of the calibrators: HD 92933 (K=5.58mag) and BD- 
03 38 26 (K=6.70mag) . We used the temperature reddening 
from lTokunaeal (|2000[ ) to convert to L', and scaled the mag- 
nitude with respect to the relative flux observed on both 
calibrator and science star: HD 92933 is 1.09 ± 0.03 times 
brighter than HD 92945, and BD-03 3826 is 1.11±0.05 times 
brighter than HD 141569. As a result, we estimated the 
L' magnitude of HD 92945 to be 5.58 mag and HD 141569 
to be 6.74 mag. We checked that these values are roughly 
compatible with the magnitude derived s olely by using the 
isochrones for pre-main sequence stars (jSiess et al.l 120001) 
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giving 5.6 and 6.3 mag, respectively. Scaling the bright- 
ness to derive the absolute magnitude at 10 pc (assuming 
distances stated in Sect. 14. 1]) . we obtained ML'=3.87mag 
(HD 92945) and ML'=1.76mag (HD 141569). 

Thus, the 5 a upper limits correspond to a nondetec- 
tion of up to a absolute magnitude of Ml' = 10.4 mag for 
HD 92945 and A/l' = 7.6 mag for HD 141569. To convert 
our detection limits to planet ary masses limits, we used 
£>t/5ry evolutionary models ()Chabrier et all |2000D con- 
volved with the NaCo filters. Towards HD 92945, the mod- 
els put a limit on the mass of a companion to 18 Mjup at a 
separation of 1.5 Au(A/-D). Towards HD 141569, the models 
limit the mass of a companion to 22 Mjup at 7au. 

5. Conclusion 

We have shown that aperture masking gives detection limits 
of the order of A L' mag — 6, with an inner working angle 
close to A/2D. These results confirm previous detection lim- 
its obtained by the same aperture masking technique on the 
Keck telescope (Kraus ct al. 2011; Hinklcv et al. 2011) In 
terms of scientific impact, this observational domain is im- 
portant because it corresponds to a few astronomical units 
at a hundred parsec, the distance where the closest forma- 
tion regions are. The scientific importance of this parame- 
ter space i s highlighted by T C ha b detected by the same 
technique (JHuelamo et al.ll201ll) . 

Simulations of the older debris disk tend to show that 
a magnitude or two in dynamic range is still needed to 
observe disk shaping planet s. Considering HD 141569, for 
example, iReche et al.l (|2009[ ) show that the disk geometry 
could be best modeled by a flyby star and a planet of a 
few Jupiter masses. For these kinds of objects, a detection 
limit of two Jupiter masses would require a precision on the 
closure phases of 0.01 degree, something only possible if we 
understand how to precisely account for the systematics 
errors. 

This paper has to be considere along with othe r direct 
detect ion techniques. In the case of HD 92945, Bill er et aD 
()2007() used a spectral differential technique to derive an up- 
per limit of 10.8 mag at 0.5" in the H band (5 cr). This limit 
assumes a clear methane signature (a spectral type T8 for 
the planet) which may not be the case. This assumption is 
not required for a differ ential technique using s ky rotation: 
in their Gemini survey, iLafreniere et al.l (|2007[ ) obtain the 
same dynamic range (10.8 mag; 5 a) at the slightly larger 
separation of 0.75" . Overall, differential imaging techniques 
often give a dynamic range that is higher than SAM, but 
with larger inner working angles (typically a few A/D). 

The capacities of future extreme adaptive optics to ob- 
serve at smaller inner working angle, compared to the abil- 
ity for aperture masking to probe at higher dynamic range 
by understanding the systematics, will decide the fate of 
SAM with respect to faint companion detection. 
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